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ABSTRACT 


In a previous paper, the PPN (Parametrized Post-Newtonian) 
formalism was used to analyze relativistic influences on stel- 
lar stability in nearly all metric theories of gravity. That 
analysis omitted all "preferred- frame” terms. In this paper, 
possible preferred-frame effects on stellar stability are 
examined and no new instabilities are found. In particular, 
we show that: (i) Although terms linear in the preferred- frame 
velocity w (time-odd terms, analogous to viscosity and energy 
generation) change the shapes of the normal modes, their sym- 
metry properties prevent them from changing the characteristic 
frequencies. Thus, no new vibrational or secular instabilities 
can occur. (ii) Terms quadratic in w do not change either the 
shapes of the normal modes or the characteristic frequencies 
for radiel pulsations. Thus, they have no influence on radial 
stability. (iii) Terms quadratic in w do change both the nor- 
mal modes and the characteristic frequencies of nonradial pul- 
sations; but in the limit of a neutral mode these changes 
vanish. Hence, there is no modification of the criterion for 
convective stability, i.e., the standard Schwarzschild crite- 
rion remains valid. 
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I. INTRODUCTION AND SUMMARY 


In a recent paper [ Ni (1973); hereafter referred to as Paper I; equations 
therein will be referred to simply as I (25), I (5h), etc.], we used the PPN 
(Parametrized Post-Newtonian) formalism of Will and Nordtvedt (1972) to analyze 
stellar stability for spherically symmetric stars in a large class of metric 
theories of gravity. We found that (i) for "conservative theories of gravity," 
current solar- system experiments guarantee the existence of a dynamical rela- 
tivistic instability; (ii) for "nonconservative theories," current experiments 
do not permit any firm conclusion about whether relativistic effects will 
actually stabilize or destabilize stars; (iii) the standard Schwarzschild 
criterion for convection is valid. The basic equations of Paper I were rigor- 
ous within the PPN framework; but when analyzing and manipulating those basic 
equations, we ignored all "preferred- frame terms" (all terms containing the 
velocity w of the star relative to the mean rest frame of the universe). This 
paper completes the analysis of Paper I by analyzing the effects of the 
preferred-frame terms. 

At the time of writing Paper I, we believed that the preferred- frame 
terms would drive vibrational instabilities of spherically symmetric stars, 
and that those instabilities, when combined with astronomical observations 
on white-dwarf pulsations, might lead to experimental limits on the "preferred- 
frame parameters" Q^, and 0^. But as a result of the present analysis, 

the above beliefs turn out to be incorrect."*’ No new instabilities of any type 

The experimental limits on (= 7A-^ + A^ - 4? - *0, ^2 + ^2 " 

and 0^ (~ - 2j - 2 - £) due to stability observations of white dwarfs, 

as quoted in Ni (1972) must be dropped in view of the present analysis. 

Although preferred-frame vibrational instabilities might still occur in 
asymmetric stars, they would not likely give such tight limits on the (X 1 s. 


i 


1 



are produced by the preferred- frame terms. Hence, the stability conclusions 
of Paper I are completely general within the PPN formalism. 


II. ANALYSIS 

In Paper I, by calculating in the PPN rest frame of the star, and by 

describing the stellar pulsation by a Lagrangian displacement of the form 
i Qt 

|(x) e , we obtained the following linearized, adiabatic pulsation equa- 

2 

tions [ eqs. I (32) rewritten in slightly different format]: 


'The notations and conventions of this paper are the same as those in Paper 
I unless otherwise specified. For definitions of the various quantities 
undefined here, see Paper I also. 


b U + £ 0 + s£. fl + d ) Z - 0 . 


w 


1W 


( 1 ) 


n 


Here , af, 0 , a?. , and a?. are linear operators that are all independent of 
"o 2 iw 2 

w Q 

the frequency of pulsation £1. In particular, 

j f = - v|f 1 + (3/ - 1) U] AP + (37 - 1) PAU + ^ vl[ 1 H (37 - 1) U ]p| 

+ PV(AU) + 2p[A4VU + W(AU) + v(ad)] , (2) 

3 

tt 2 i = (| «2 + \ “3 - | V - \ a 2 w7wSp ? (AU 75 ' l 5 7S AU) ’ (3) 

w 

In Paper I, there was an error in the Ap/p term: i.e., [1 + -^(q^ + 

O 

- Q^) w + (37 - l) U] should be corrected to [ 1 + (37 - 1) U], as has 
been done above. 
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X ft 

£ i = i[| C^p (AU - i • VU) w - g P(a^ - 2o^) w VV^ 


+ « 2 P 


sf 2 £ = at + 
a ~ 

Also 


(w • V) V(AX - 6 • vx) - -g a x (w • l) PVU] > 

■|{a 2 - 5 X + 1 ) p(v - w) + (57 - 1 ) US - !(«! + **r + ! ') y 


Ap* = - P* div | 


00 

( 5 ) 

( 6 ) 


Notice that the linear preferred- frame terms are embodied in the operator 

J> and the quadratic preferred- frame terms are embodied in a?. 2 - [Recall 
iw’ w 

(Paper i) that w is the velocity of the star relative to the mean-rest frame 
of the universe. ] 

In Paper I, we derived stability criteria by analyzing the simplified 

pulsation equation U + il 2 =£ J i = 0. Here we look into possible modifi- 

° Or ~ 

cations of those criteria by the presence of =£ iw and =£ 2 - Terms m =£ lw have 

w 

imaginary coefficients; therefore they (like viscosity, energy generation, 

and radiative transport) might possibly affect the vibrational and secular 

stabilities of the star. Terms in =£ 2 have real coefficients; therefore 

w 

they might possibly affect the dynamical and convective stability of the 
star. But, in fact, no such effects occur. To see this, we proceed as 

follows. 

Equation (l), when supplemented by the boundary conditions and the 
expressions for the various Lagrangian changes in terms of constitutes 
a characteristic value problem for fl. In Paper I we have shown that (i) 
this characteristic value problem is self-adjoint if, and only if, 

(?) 


*1 - 2 ~ ^ ’ £2 *" £3 ^ h " ^ y 


a, = a. 
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but (ii) in the general case with w- terms deleted, although the characteris- 


tic value problem is not self-adjoint, a variational integral can still be 
constructed. 

When the w- terms are included, one can use the same argument as in 

Paper I [the passage between I (43) and I (53)] to show the following: 

* 

Take the full post-Newtonian characteristic equation (1); dot | into it; 
and integrate over the star. The result 

fi 2 J* I* • 2 £d 3 x + Q J I* • ^. w |d 3 x + ]*!*• U 0 + =? 2 ) |d 3 x = 0 , (8) 

n ~ ~ ~ w 

is a variational principle for the post-Newtonian normal modes. 

The changes in squared frequency due to preferred-frame terms can be 
obtained, to post-Newtonian accuracy, by inserting the Newtonian character- 
istic functions into the variational principle (8) and solving for ft: 

r t* ■ U o + te. ) |d 3 x 

p ~ w 2 LW ~ 

(fcft ) £ A * = - 5 — i — — — * (9) 

preferred frame r* ,„i2 ,3 ' 

J p[!| d x 

The change 3ft due to is proportional to 

1 I ' = P iw! d3x ' (10) 


Being linear in the constant vector w, raises and lowers the spherical- 

harmonic index of a characteristic function by 1 — thereby producing a new 
function orthogonal to the original one. (For proof, please see the Appen- 
dix. ) : 


3ft 


J ( { 

J (‘ 


mode ) 
mode) 


• (l - mode) d 3 x 
iw 

■ j a[ ( l - 1 ) - mode ] + b[ ( £ -t 1 ) - mode ] 


d°x 


0 -(H) 
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Thus the 3 ? term has no effect whatsoever on the pulsation frequencies of 
’ ‘ iw 

the star to post-Newtonian accuracy. Hence, to the same accuracy the stabi- 
lity criteria are not affected either — and, in particular, there is no 
post-Newtonian vibrational or secular instability. 

Might the preferred- frame velocity w drive a vibrational or secular 
instability at higher orders in the post-Newtonian expansion (post-post- 
Newtonian, etc.)? There are two ways that it might do so: (i) by higher- 

order corrections to the imaginary operator 3 ^, and (ii) by higher- order 

effects of the original *?„ (eq. [4]). Although one is not now in a position 

iw 

to evaluate the first way,, it is clear that the second way cannot produce 

an instability: Use higher-order perturbation theory to solve equation (l) 

to higher-order accuracy in a?.. . Since / is purely imaginary and af, Q , =£ 2 , 

_p are a ix real, every term in a perturbation series expansion for Im(fl) 

must contain ^ an odd number of times. [This can be proved by an actual 
"iw 

expansion calculation despite the nonself-adjoint nature of the problem.] 

But iP only raises and lowers the spherical-harmonic index of a character- 
“1W 

istic function by 1 , while £ , £ , and only raise and lower the spherical- 

w ft 

harmonic index by 2 or do not change it [for proof, see the Appendix]. Hence 
every term must contain an integration of the form 

f ( t - mode)* (£' - mode) d’x , * ( 12 ) 

where one of l and V is even while the other is odd. By orthogonality of 
the different spherical-harmonic modes, ( 12 ) vanishes and 

Im{Q) = 0 . (13) 

Therefore no higher-order vibrational or secular instabilities can occur. 

1 
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By equation (9), the change in squared frequency due to the rela- 
tivistic j£.p operator is 

w * ■* 

5il « i I 2 6 d x • W 

W 

The first term in equation (3) for sf ^ | can be combined with the term 

w 

5V(AU) in s£ q I [eq. (2)] to give a renormalization of the gravitational 
constant 

G -> [1 + a 2 + -| a 3 - ~ a x ) w 2 ] G . 


Since Newtonian stability is not sensitive to the value of the gravitational 

L 

constant, this term does not affect stability (even though it does affect 

the values of nonzero pulsation frequencies). The second term in J ^ | 

w 


[ eq. (3)] can be put into the form 


where 


and 


\ “ 2 » r » 5 PV(5U 7S - \ 8 8U) , 


5U^(x) = J 5p(x ? ) 


(x - x' )(x - x’ ) 

\ r-tr r%M t \ O O' 


7 7 5 


x - x 


. .3 


dx' 


(15) 


( 16 ) 


5U(x) = J 5p(x' ) 



(17) 


The symbol 3 denotes the Eulerian change in the quantity that it qualifies. 
For radial pulsations 


8U 7B (x) = | 8 75 & u (x) ; (10) 

therefore equation (15) vanishes and there are no changes either in the shape of the 
normal modes or the characteristic frequencies. For nonradial pulsations, 
although the normal modes and characteristic frequencies are modified in 
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general [see the Appendix for an explicit expression for 5 ft ], at the onset 
of convection 

8 p(x') = 0 , (19) 

so 5 U = gU = 0 and equation (15) vanishes. Hence by equation (lb), neither 

75 

spherical stability nor the onset of convection are affected by a? Q . 

w*~ 

III. CONCLUSIONS 

Preferred- frame terms have been shown to not affect stability criteria 
for stellar pulsations about spherical equilibrium, even though they do 
modify the characteristic functions and (in some cases) characteristic fre- 
quencies. This enables us to conclude that the stability criteria found in 
Paper I are completely general and can be used for all metric theories of 
gravity which fit into the PPN formalism in their post-Newtonian limits. 

I am much obliged to Dr. Kip S. Thorne for many helpful suggestions 
and critical comments. 
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APPENDIX 


SOME DETAILS OF THE CALCULATION 

In this appendix we will derive an explicit version of equation (9) 

for and will prove the statement ("S”) that only raises and lowers 

the spherical-harmonic index of a characteristic function by 1, while =£ 2 , 

w 

f and af. ^ only raise and lower the spherical-harmonic index by 2 or do 

° n 

not change it. 

To accomplish this, we analyze the Lagrangian displacement £ into 
normal modes belonging to different vector spherical harmonics, i.e., 




= 


■ 1 dx(r) 

j l(t + 1) r dr 


a© 


(Al) 


(A2) 


and 


1 

- l{i + 1) r sin 0 


dx(r) ,<■»(©,*) 
dr 00 


(A3) 


(;: ^ are physical components, not covariant components). For con- 

-^gj^iaace, we choose our coordinate system such that the z— axis is in the 
direction of the "preferred-frame” velocity, i.e., w - wz where z is the 
unit vector along the z-axis. 

2 

Without the preferred- frame terms and s£ 2 , the operator s£ q + H ■£ , 

w w H 

is invariant under rotation; hence the characteristic value problem has m- 


degeneracy in this approximation. With the preferred- frame terms, the 
degeneracy in m is broken; but since the problem still possesses azimuth 
symmetry, normal modes have definite m-numbers. Therefore (Al)-(A/, ) are 
normal modes for the characteristic problem to the first approximation, and 
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can be used to evaluate the changes in characteristic frequencies due to 
preferred-frame terms by equation (9). 

To evaluate equation (9) and to find the effects of and ^ on 

iw ' 2 

w 

normal modes, we first derive a number of formulas involving spherical 
harmonics. From the integration formula involving three spherical her- 
monics [cf. e.g, , Edmonds (i960) p. 63] 


J' Y Vl (9 ’^ d « 


(2£ 1 + l ^ Zt 2 + l >< 2f 3 + 1 >l l/ Vl *2 A *2 * 


^0 0 0 / \1 m 2 m 3j 


, (A4) 


we have 


J* Y Y Zm ( ' 0 ’V ) COS ° dQ 


(as) 


and 


S Y ae 


sin 0 df2 


- (- l) m [2(f - m)(£ + m + l)(2 V + l)(2f + l)] 1 / 2 


(V l l\ / VI 


^2 ^3\ 

where I ) is the Wigner 3-j symbol. I Q ” ~J is nonzero only if 


0 0 0/ y- m 1 m + 1 - 1, 

V St l\ 

igneL o-j symDOi. i 

1 = f + 1 or I' = f - 1, 

From the gradient formula [e.g., Edmonds (i960) p. qU] 

i\ 1/2 I 


V[*(r) Y M (e )t p)] = - 

( t x 1 / 2 

+ \2f + 1 / 


A * + 

dr H ' r 


(/r - 

~ J $(r) 


f £+ lm 


A cp) 


Vl J°’V> ’ 


(A7) 


, (A6) 
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and from the definition of vector spherical harmonics, we derive 
z • V [®(r) Y ^ m ( e ^)l = 


( £ + m + l) ( £ - m + l) 

1/2 

/jL A ' 

(2£ + 1)(2£ + 3) 

- _ 


\dr r 


£-i-lm 


' i 

(£ + m) ( £ - m) 

1/2 


( 2 £ + l) ( 2 £ - l) 


(dr + r f 


£- lm 


and 


J-jS • v (*( r ) Y * £m ) 
( £ + 


z * V 


Y iJ! 

m !!!;.-.•=; i; [& - 7] [(£ - 7 
[fe * ht] [te * ^y 


•(£ + m) (£ - tn) |/d ^ + 

+ C 2/ t 'f) 7 27 vry + r 


From the formulas 


d/r^ + r |g - 2 rr* cos' 7 

c^r 


r - r 1 cos y 


[2 ,2 

/r + r f 


2 rr f cos y 


and 


/ 


2 2 

r + r' - 2 rr' cos y 


OO 

^ £h-1 a £ 

!=0 r’ 


we find 


/ 


2 2 00 
r + r T - 2 rr’ cos y - E 

£=0 


£+2 


1 


2£ + 3 , £+1 2£ - 1 , £+1 

r r 


(A 8 ) 


(A9) 


(A10) 


T7T p „( cos ?), ( r < r ') (ah) 


p « (cos y). (A12) 


By the addition theorem, we can then write equation (A12) as 

£+2 


lx - x' | = 4jc E Z p j 
£-0 m--£ 


- 

1 1 

+ 12 ! + 


3 , £ + 1 

r 


- 2tH jrrf'J 0 '’*') Y zJ 0 ’^ • 

Now by (Al)-(A13) and by the expressions [cf. I (38), (39), (1+0)] 


(A 13) 
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&U = - J div[p(x f ) e(x' )] 1 dx' , 


X - X’ 


(AUO 


&U 75 = " J div[p(x' ) $(x' )} 


(x -x ? )(x - x ' ) 

x 7 7 5 & 

I . ! 3 

x - xM 


dx’ 


(A15) 


and 


5X = - J div[p(x' ) |(x' )] |x - x’ | dx' , 


5U 7B = 8 7B BU ' ’ 

76 70 9x 7 Sx u 


we can readily verify, by straightforward calculations, that 

(i - mode) = a[(£ - l) - mode] + .b[ (f -i l) - mode] 


(Ale) 


(A17) 


(A1Q) 


and 


(j _ mode) = c[ (i - 2) - mode] + d[ f - mode] + e[ (l + 2) - mode] . (A19) 


w 


Equations (A18) and (A19) together with the invariance of and =£ 2 under 

rotation prove the statement "S" made at the beginning of this appendix. 

2 

To derive an explicit formula for 5 ft [eq. (9)], we note that 

J“ 5* A = 0 

[by (A18) ] . Therefore the integral remaining to be evaluated is 

J* e* • * 2 1 d 5 x. . 

w 

By (Al)-(A17) and the formula [by (Al)-(A3)] 


V • (p|) = ”2 




(A20) 


we have, after some calculations and reductions, 
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X E* • £ 2^ d 


2 * „ 2 
QLw 


• 2 ^“ X “ 21+1 ~2 
w 


JT 


d(pt) 


(l + m 4 - l)(l - m + l) /d 2 2 

(2 1 + l)(2i + 3) \ dr 2 + r 


dr’ 


dr 


- o Jfc-1 

fd(p^) 


p dr J 

[dr " M dr 

2(2 + 2)\ 

2+2 

/ 1 ^ 

r 2 ' 

\22 + 3 

2+1 




I 


> 


21 


r 1 v 

ViTirr)- 


( i + m) ( 2 - m) / d _2 d 

(22 + l )(22 - 1 ) 


> 


dr 


2 r dr 


(j - i)(* + 


1 + 2 I I 

!< . 

• /\22 + 3 2+1 2.2 -1 j2- 1/ 


dr dr ' 


(A21) 


> 


> 


In general (A21) does not vanish. But for 2=0, (A2l) vanishes; and for 
1 ^ 0, at the onset of convection where, $p = 0, i.e., where 


d(p^) 

d r 


P 


dx 

dr 


0 


(A22) 


(A2l) also vanishes, in agreement with the results in §11. 
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